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The influence of spin on a photon's motion in a Schwartzchild and FRW spacetimes is 
studied. The first order correction to the geodesic motion is found. It is shown that 
unlike the world-lines of spinless particles, the photons world-lines do not lie in a plane. 

1. Introduction 

The motion of spinning particles, test ones or tops, in a curved spacetime has 
been a subject of great interest for the last nearly 80 years. It is amuzing to note 
that, even today, there is no general consensus on the behaviour of particles with 
spin in external gravitational fields. It seems that the first to discuss the motion of 
free extended material particles with spin in an external electromagnetic field was 
Frenkel 1 . The classical test body was considered to be so small compared with the 
background curvature length-scale that all its multipoles beyond the dipole could be 
neglected. The equations of the free motion of such a 'pole-dipole' particle for the 
linearized theory of gravitation were obtained by Mathisson 2 . Papapetrou 3 worked 
with non-singular localized energy-momentum tensor Ty for a spinning particle (in 
the limit of a vanishing mass) and obtained equations for moments of T^ . Dixon 4 
found a way to make Papapetrou's type of argument covariant at each step of the 
derivation and also found equations which referred to paths described by an arbi- 
trary parameter rather than the world-line length. He came up with the equations 
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of motion in the form 



Dr> 1 1 

Dq~2 K * km ' 
DS 11 



, (1) 

with q an arbitrary parameter along the curve, v l = dx l {q)/dq, D/Dq a covariant 
derivative along the world-line x 1 (q) and p % and S tJ defined as in Papapetrou 3 . The 
explicit form of generalized momentum p, containing the additional term depending 
only on spin and connection, was shown in [5], in complete agreement with the 
definition of Papapetrou. In the derivation of these equations by Dixon nothing 
required the object to have a rest mass or required g i jdx l dx : > to be nonzero. They 
were, therefore, sufficiently general to relate to any localized object. 

However there have not yet been satisfactory direct derivation of these equations 
for massless particles with spin. In [15] Baylin and Ragusa treated massless parti- 
cles as traceless, but the authors couldn't come out with the final answer whether 
"traceless" -massless particles' motion deviate from null geodesies. At the same time 
they mention that "we have been unable to prove that null geodesies necessarily 
follow for the general case of spin in curved spacetime. It even seems unlikely in 
view of complexity of [Papapetrou-Dixon] equations" . 

Explicit coupling to curvature actually turns out to be the general rule in the 
field theoretical case, so that strictly geodesic behaviour is not to ever be expected 
in general. Because fields are defined at all spacetime points and not just along the 
trajectories associated with point particles, the associated second order field equa- 
tions are sensitive to the values of the connections at differing points. Both Dirac 7 
and Maxwell 8 ' 9,7 equations lead to an implicit expressly non-inertial coupling to the 
Riemann curvature tensor. And even though it is possible to remove the Christoffcl 
symbols at any given point (or even along a single curve), it is nonetheless impos- 
sible to remove them from an entire region, with the curvature dependent terms 
we find in the second order wave equations being the field theoretic generalization 
of the (covariantly describable) geodesic deviation found for pairs of nearby falling 
particles. 

In the previous paper (paper I 10 ) we have derived the equations of motion for a 
vector particle with spin from a Lagrangian. We have worked out the approach that 
gives a satisfactory approximation to the wave as some curves that can be called 
"rays" and, at the same time, takes into account the polarization. The obtained 
equation coincided with the Papapetrou equation and is independent of a mass of 
a particle. In the present work we apply equations derived in paperl for photons 
transversely propagating in a Schwartzchild metric. We develope the technique to 
calculate the corrections to the geodesic motion. We show the deviation from a 
geodesic in Schwartzchild metric and examine the consequences of the deviation in 
different gravitational backgrounds. We present a general discussion of the results. 



2. Geodesic congruences for Schwartzchild metric 
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2.1. Hamilton- J acobi Equation for photons 

The Hamilton- Jacobi equation (HJE) for photons reduces to the eikonal equa- 
tion, 

dip dtp 



In Schwartzchild geometry 



2m\ 2 / 2r7j\ ' 



ds z = ^1 - dt l -\\- '-J dr A - r l d6 z - r l sin^ dip' , 

with the coordinates: x^ — {t, r, 8, ip} = {0, 1, 2, 3}, the metric coefficients are: 



^ = diag(-i— , -(1 + 2^,-1, -^-^-) , (+) 
1 + 20 r z r z sin r 



where by <p we denoted the "Newtonian" potential <f> = —m/r. We employ here the 
Planck units c = h = G = 1. Eq. 2 becomes 



1+2 nf) 2 - (1+2 «(^) 2 ^(Sr^(l) 2 -- <4 > 

Following the usual approach to separating the variables in HJE, we represent the 
desired solution in the form 

il> = Et + L x <p + R(r) + Q(6) , 

with the convention that \fa? = ±a. Substituting it in Eq. 4, we obtain the equa- 
tions where variables have got separated indeed and produced an arbitrary constant 
(the constant of separation, total angular momentum) L: 

(1 + 20> 2 [R'f = (e'f + J i~ = L\ 
1 + 20 sin 

from where we find functions 6 and R, and get the solution in the form 



sin 2 i 



2.2. Associated Tetrad Basis 

Any coordinate system {x 1 } in the spacetime specifies the natural vector {di} 
and covector {dx 1 } frames at each point. We defined the metric by the scalar 
products of natural covectors 



g' 3 = (dx J ,efaJ) , 
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because our main requirement to the metric is that it shall separate the HJE for 
geodesies, which looks much more simple in the contravariant components than in 
covariant ones. 

We choose for simplicity the equatorial plane, 9 — tt/2 (L = L z = p v ). Conse- 
quently, the dynamic phase ip becomes 



f JE 2 r 2 - L 2 (l + 2(f)) 

where we fixed the signs as 

tp + = Et - R + Lip , 
ip- = Et + R — Lip , 

and write the gradients of the phase 

#± = Edt ± Ldp =F R'dr . 

Now we can construct the associated basis on which the connection form is identi- 
cally zero. Taking the scalar product {dtp +, dtp-}, we obtain 



(#+,#_) - E*(dt,dt) - L z {dip, dip) - {R'f (dr,dr) = , 

since from the metric coefficients (Eq. 3) we have (dt,dt) — g u = 1/(1 + 2(f)) 
(d9,d6) = g 66 = -r- 2 , (dr,dr) = g rr = -(1 + 20) and (dip, dip) = g w = 
— 1/r 2 sin 2 9. The first pair of basis vectors is thus 



2E 2 



dip± #±\/T + 2(f) 



y/(dip+,dll)-) EV2 



Since is a null vector, we need to add timelike vectors to our basis. Applying the 
Hamilton- Jacobi theorem 11 to the solution (Eq. 5) we can construct the integral of 
motion in the form of a function whose values are constant on the trajectory 

8L V 9L 

We choose 

dV = dp — dr . 

oL 



The scalar product (dV, dV) is 



(dV,dV) = (dip,dip) + ( — \ (dr,dr) 



dL J " ' ' r 2 - D 2 (l + 2(f)) ' 

where we denoted L/E = D, with D having the meaning of the impact parameter. 
Thus, 

v 2 = ■ dV = dVJr 2 -D 2 {l + 2(f>) . 
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As a fourth basis covector we choose d9: 



v 3 = 



= rd9 . 



Thus, we constructed the associated basis of 1-forms, 



v ± = [Edt ± Ldip t R'dr] 



v 2 = ^r 2 - D 2 (l + 2<j>) 
u 3 = rdO . 



A dR ' A ' 
dtp —dr 



(6) 



This basis is orthonormal in the sense of (^ ± , ) = and (v*,^) = 1. The basis 
dual to {v a } is the vector basis {n a } (v a (n}>) = Sg) 



n± 



n 2 = 



n 3 



E d 9 L d 

J — W) - T R(l + 2<j>)-±-,- 



EV2 

D(l + 2<p) d , y/r 2 - D 2 (l + 20) d 
r 2 dip 



dr 



+ 



(7) 



1 d 



r 89 ■ 

We also introduce orthonormal sub-frames, constructed from the associated basis, 



v 



.0 



v 



and 



n 



\/2(T+20) 9t ' 



EV2 



L _d_ 

dip 



£-{l + 2<t>)R>£ r 



According to the definition Def. 1 
tions of the vector basis are 



h L = 1 

V 2 ( 1+2 ^) 



{n h } : n b = hi 



the non-zero tetrad func- 



h r _ -P(l+20) 



K = -^^r>-D\l + 24>) 



'1+20 



hi 



2.3. Polarization basis on the geodesic 
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We assume a narrow beam of rays and neglect transversal derivatives — in the 
zcroth approximation polarization basis is propagating parallely along the ray. We 



can write tensors as F— F ab f ao , where {f a b} — tensor basis 



fab = n [a <g> n b] = i (h l a h{ - h l b hij 



d d 



(8) 



and square brackets mean antisymmetrization. We define vertical (v) and horizontal 
(h) linearly polarized waves 

A v = A 2 ri 2 = (a simp) ri 2 , 

A h = A 3 n 3 = (acostp)n 3 , (9) 

where a is an amplitude. The waves with left and right circular polarizations (LCP 
and RCP, correspondingly) are their linear combinations, 

= A v + Ah = (a sin ip) n 2 + (a cos ip) ri 3 , 

= A v — Ah = (a simp) ri 2 — (acosip) n 3 . (10) 

The derivatives of (9) are 

? , ; m - ( aE ^ \ - 

A v = A n 2 = (atp cos ip)n 2 = —^=^= cos ip \ n 2 , 

V VI + 2-7T / 

A h = A 3 n 3 = -(aipsinip)ri 3 = - ^ ^^^ smip^j n 3 , (11) 

where we used the fact that on the geodesic the connections are zero and basis 
vectors ri a are parallely transported (by construction). The operator 'dot' was 
defined in the following way. By the definition, a gradient of a function in dual 
frames is 

df = (n a o /) v a . 
We take the derivative only along the beam, hence, 

df = (n_ o /) v+ , 

and replacing v + from Eqs. 6, 



At t\ V 1 + 2^ , ,+ 

df = (n_ o /) ——j=-d%p + . 



Finally, the 'dot' operator takes the form 
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The spin tensor is denned as (see Appendix A) 



?ab = \j 



A 2 A 3 - A 3 A 2 



d 3 x . (12) 



The only non-zero component of the spin current 

a 2 E 



J 23 



^2(1 + 20) 



and the only non-zero component of spin for the LCP polarization we find from 
Eqs. 10, 11 and 12 

L »S'23 = / , a = = (cos 2 ib + sin 2 ib) d 3 x = [ , = = = d 3 x = const = 1 . 
J ^2(1 + 20) V ' J V2(l + 2^) 

For RCP, correspondingly the sign is opposite, 

R 5*23 = — 1 • 

It is convenient to introduce the following notation for the spin-gravity coupling 
term, the tensor 

Tji ^ pi Qkl 

n o - ^ m ' 

where in order to express the spin tensor in the coordinate frame, we use the tensor 
basis defined in Eq. 8, 



5 23 - ±1 = S 23 = S U h%hf = ± (/, 

and 

s kl = S 23 (f 23 f . 

From six components of {/23} basis only two survive 

(723 j - -2 h 3 f h ^2 ' 



23 

kl) 



i/23) - -2 h 3 h 2 ^3 ' 

where the tetrad functions are found from Def. 1). 
Thus, for LCP the components of this tensor are 

L tt9 1 r>9 (t xflr _ mD 

H r - 2 H r9r(h3) - , 



L tt8 1 r>8 (f ^_ m v /r 2 -/J 2 (l + 20) 

^ V -2 K ^U23j - 7^1 
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where the corresponding components of the Riemann tensor arc 

r>0 m 

-K rftr = ~ 



r 3 (l + 2</>)' 

U V8<P r ■ 

3. The Papapetrou Equation 

The first equation of (Eq. 1) can be written in the form 

V (13) 

With our definition for the circular polarization the equation transforms to 

p i = ±H i j x i (S 23 ) (14) 

with '+' corresponding to LCP and '-' - to RCP. The geodesies we are considering 
lie in the equatorial orbit 9 — ir/2 and, hence, 9 — 0. Since the only surviving 
H tensors contain 9 components, in the Eq. 14 wc have non-zero right hand side 
only for ^-component of the momentum. We are considering the deviation 8x l (i})) 
from the geodesic x l (tp) as our first approximation, Sx l (ip) = (0, 0, 56(ip), 0). The 
Papapetrou equation for all coordinates but 9 coincides with the geodesic equation, 
and for the 9 coordinate we have: 

P 6 = 0, 
Sp e = ±H e jX j , 

where in the zeroth approximation (geodesic motion) the spin-dependent term of 
the generalized momentum p l vanishes. To find the components of the velocities 
we need to find the geodesies in the explicit parametrized form. If we express the 
solution as it>{x l , a a ), where x 1 represents variables and a a represents constants, 
then equations 

dip 

- — = ip a {x\ a a ) = const . 
da a 

give the analytical representation of geodesies in parametrized form (Hamilton- 
Jacobi theorem 11 ). Applying this to the Eq. 5 gives geodesies in the form of the 
following system of equations with r as an independent parameter: 



dE 



to ■ 



dtp± 



t n = ± f Erdr = t(r) ■ 

J (l+20)^/£ 2 r 2 -L 2 (l+20) V ; ' 



W - U> = ± f — ML = tp(r) . 

* YU J ry/E 2 r 2 -L 2 ( 1+24,) ^ V ' 



(15) 
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Now we need to find the relation between momentum and velocities. The compo- 
nents of the momentum vector are obtained from the solution of HJE, 



Pt = = E ; p r = R' ; p v = L ; p e = , 



and, accordingly, 



p* = -A_; p r = -(l + 2 l p)R'; p* = ~; P 6 =0. 
1 + 2tp r l 

Now we find (x l )'-s from (15) 

(**)' = Er = ; (x r )' = 1 ; (x 9 )' = ; 

(l + 2^) v /£ 2 r 2 -L 2 (l + 20) 

(x v y = —= L = , (i6) 

where 'prime' denotes differentiation with respect to (w.r.t) r. By taking out the 
common multiplier, we express the momentum through velocities 



pi = v /g2 r 2_ L 2 (1 + 20) ^ = ^-^(1 + 20) ^y ^ 

In our first approximation we can assume the same relation between the velocity 
and the momentum as on the geodesic, 



5p e = E ^-nM + 2* Uel 

r 

By 'prime' we mean the covariant differentiation with respect to r. However, we 
can express Sp e through our associated basis, 



Sp e = ^/~^8p 3 = -Sp 3 . 



(18) 



According to the zero-connection property of the associated basis, we can replace 
the covariant differentiation of Sp 3 w.r.t. r with ordinary differentiation w.r.t. r. 
The same would hold for a product of Sp 3 with any scalar function, as in (18). 
Hence, we write the Papapetrou equation as 

The contraction with the velocity gives the right-hand side of the equation 
with (x J )'-s from Eq. 16. Thus, the Papapetrou equation is 
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and, according to (17), 



Ey/r*-D*(l + 2<l>) \ _ mD 
r I 4r 5 



Thus, the trajectory of a photon world-line is: t(r), ip(r), 6(r) = ir/2 + 69(r), where 
for LCP (+) and RCP (-), correspondingly, 



±S9(r) 



iD 



16E 



dr 



r 3 ^/r 2 - D 2 (l + 2(f>) 



(20) 



with gravitational deflection term ip(r) defined in (15). With (1 + 2<fi) = (1 — r g /r) 
and m = 2r g , Eq. 20 takes the form 



dr 



Vr 2 -D 2 {l~r g /r) 



(21) 



This formula shows the effect of dispersion, the energy-dependent deviation from 
a geodesic motion due to the spin. To see whether this quantity will have any 
appreciable effect on the propagation of light, we have to make an estimate of the 
order of magnitude of this integral. It is more convenient to express the result in 
terms of the distance ro of closest approach, rather than the impact parameter D. 
At r 



r(l + 2<p) 



r 



0. 



Since the difference between D and ro is small, we write D = tq + 5ro, and inserting 
it into the previous equation gives to the first order <5ro w r 9 ro/2(ro — r g ) and 
D w (2r + r g )/2. Using ro instead of D in the integral (21) gives 



±S6(r) = - 



8E 



r r g dr 



'8eJ 



r r g dr 



1 + 



(2ro+r a ) 2 
4 

r a r g 



r 3 v^^V 2r(r 2 -r 2 ) 



Up to the first order in r g /r we obtain (E = 1/A) 



±S9(r) = 



r_g\ 
I6r 2 



ro^r 2 - rl 



And at the observer the deviation from the 9 = tt/2 plane is 



£50(+oo) = T 



16r 2 



Motion of a vector particle ... 11 



We notice that at r g — ► 0, or at E —* oo, the spin-dependent deviation vanishes. 
Spinless particles always move on plane orbits, in the case of a circularly polarized 
photon with the spin (helicity) parallel to the velocity, the perturbed motion is not 
plane. 

We feel that we shall comment here on the problem of the deflection of light. 
Contrary to the reports 17,18 about the polarization dependence of light deflection in 
a Schwartzchild metric, and in concordance with the pioneering work of Corinaldesi 
and Papapetrou 16 , we conclude that spin gives no contribution to the deflection (the 
deflection (Eqs. 15) is identical with that which is found from using the equations 
of geodesies). 

Of course, the result with X/tq has an extremely small value and it is evident 
that the effect certainly has no observational consequences in the gravitational fields 
appropriate to present-day astrophysics and cosmology. Nevertheless, we feel that 
the main significance of the result is theoretical and does not depend on its ob- 
servability. Moreover, in view of recent findings of the appearence of a chaotic 
behaviour in the motion of a spinning particle in a Schwartzchild spacetime due to 
the spin-orbit coupling 13 , and of a report of a measurement of a time delay between 
RCP and LCP signals from a pulsar PSR 1937+21 19 , we believe that a study of 
spin effects on the orbital evolution of relativistic systems is very important from 
the viewpoint of observations as well as from an academic one. 

3.1. Robertson- Walker Background and Radial Motion 

From the general form of the Eq. 13 it is clear that in the flat spacetime (R^ Kp = 
0) or for the spinless particle (S^ = 0), the equation reduces to the geodesic 
equation or, in case of a massless particle, to the null geodesic equation. 

In general metric at every point of a spacetime there exists a Weyl principal 
tetrad in which the Riemann tensor takes its normal form 12 (it is strictly diagonal 
in pairs of indices). If the direction of the momentum of a spinning particle coincides 
with one of the principal tetrad directions, then due to the Tulczijew constraint 14 
(one of the supplementary conditions to the equations of motion, introduced in 
order to reduce the number of spin independent degrees of freedom, p a S a ^ = 0) the 
Papapetrou force is identically zero. In the homogeneous gravitational background, 
or background with constant curvature, any direction is the principal direction, and 
hence the Papapetrou force is always zero in such metrics. Since the FRW spacetime 
is isotropic and homogeneous, the origin of polar spherical coordinate system, for 
example, can be any arbitrary point. If we place the origin on a geodesic, the 
geodesic becomes strictly radial in this coordinate system. Thus, the photon's 
motion is geodesic in FRW spacetime. The same can be shown for any spacetime of 
constant curvature, for example, de-Sitter, the fact first discovered by Bailyn and 
Ragusa [6]. 

The radial direction in Schwartzchild metric is both principal (in the sense of 
Weyl principal tetrad) and geodesic and, thus, there is no deviation. The motion is 
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just pure radial geodesic with 



f dr 
t(r) =t ± J - 



+ 20 ' 
Edr 



0(r) = Const , 
ip(r) = Const . 
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Appendix A Spin tensor derivation 

In paper I we have defined the spin tensor as following 



mS ab = i (A a A b - AbAa^j . 



However, it was pointed to us recently" that in that definition the dimensions are 
not correct. Removing the mass term restores the correct dimensionality. It should 
be noted, however, that neither the derivations of the equations of motion for spin 
and momentum nor the final results of the paper I suffer from this redefinition. The 
derivation of the equation of motion of spin does not depend on the mass term, and 
in the last equation for the momentum (Eq. 31, paper I) the mass term is absorbed 
in the momentum term 

9ijm ~Di = 9l] "757 = x jil k • 

We will present here the proper derivation of the spin tensor which was just pos- 
tulated in paper I. Using the Noether theorem we will derive the spin tensor as a 
current conserved under the rotation of a local orthonormal frame as an operation 
of internal symmetry (see, for example, [20]). The field Lagrangian is 

2£ F = -A 2 + m 2 A 2 , 

where 'dot' means covariant derivative. Under the infinitesimal rotation of the 
frame 5a%, the new Lagrangian is the function of 

Cf = £f (A, A, da, See) , 

where 

5 A = 5a ab S ab A . 

a We express our thanks to Dr. M. M. Sheikh-Jabbari from IPM, Tehran, Iran. 
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Thus, the new Lagrangian is 

2£ F = C F + 5 (-A 2 + to 2 A 2 ) 

Let us compute it, 

C F = -l(A + 5A)- 2 + ^(A + 5A) 2 = 



£ F - A 
The variation is 



-\ (A + (SA)') (A + (5A)-) + (A 2 + 2 (5A) A) 
(<Sa ab S Qb A)' . 



<5£ F =- (As a6 A) (<Sa ab )* . 
The variation of the action integral 

= 61 = J d 4 xSC F = - j d A x (A£ o6 A) (Sa ab )' . 

Evaluating this integral by parts and assuming that at infinity the fields go to zero, 
the variation is 

51 = J d A x (Sa ab ) (AS^A) " , 
Using Noether theorem 20 , we define the current 

A5 afc A = (j ab ) . 



In the world-tube, where we defined the field, the spin tensor will be 

S ab = J (jab) d 3 X . (A.f) 

Since the spin of the vector field is I by definition, we do not need to integrate (A.l) 
over the whole space and can express spin directly through the orthonormal basis, 
defined along the world-tube. Thus, since our spin matrix is 

{Sab)cd = J (VacVbd - rjadVbc) , 

we obtain the expression for spin current as 

A c A d (S ab ) cd = l - (A a A b - A b A a ) . 

In the world-tube we assume only plane waves and, thus, identify the spin current 
with the spin. Thus, our expression for the conserved quantity, spin of the field, is 



S ab = \j (A a A b - A b A a ^j d 3 x . 
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We shall note that the spin tensor S a b is defined now only on a geodesic, in the same 
way as velocity and momentum vectors are defined, which completes our transition 
from field description to particles. 
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